A FOUR-VERTEX, QUADRATIC, SPANNING FOREST POLYNOMIAL 

IDENTITY 
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Abstract. The classical Dodgson identity can be interpreted as a quadratic identity of 
spanning forest polynomials, where the spanning forests used in each polynomial are de- 
I fined by how three marked vertices are divided among the component trees. We prove an 

t-H analogous result with four marked vertices. 

o 

o 

Q 1. Introduction 

Let G be a graph with m vertices, n edges and let the i th edge be assigned a variable a,, 
i — i Then we define the graph polynomial of G as 

fH TCG e$T 

where the sum runs over the spanning trees T of G. The reason why we pick edges not in 

the trees is that this form arises naturally in quantum field theory, see for example (3J El [10] • 

We can also obtain this polynomial via the matrix-tree theorem. Let A be the nxn diagonal 

matrix with the variables aij. Orient the edges in the graph and let E be the signed m x n 

incidence matrix for this orientation. Let E be the matrix E with any row removed. Define 

^sO the m + n by m + n block matrix 

00 

M G A ET1 



(N 

o 



^ -E 

Then the matrix-tree theorem states that 

m G = det(M G ) 



To put this in the usual form of the matrix-tree theorem, note that A is invertible, so we 
can calculate the determinant using the Schur complement; in this case, 

det(M) = ax ■ ■ ■ a n det(0 - (—EA~ 1 E T )) 
= a x ---a n det{EA- l E T ) 

where EA~ l E T is the graph Laplacian with a row and column removed and with inverted 
variables. See also Proposition 21 of [1]. 

There are two important ways to generalize - one via the polynomials and one via the 
matrix determinant. Let P = Pi U • • • U be a set partition of a set of vertices in G. Then 
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define the spanning forest polynomial for G and P as 



FCGe^F 

where the sum runs over spanning forests F of G composed of tree components T 1; . . . , 
where the vertices Pi are in tree Tj. Alternatively, let J, J, be sets of indices with \I\ = | J|. 
Define the Dodgson polynomial as 

^ = det(M G (J, J)) x 

where Mg{I, J) is the submatrix obtained by removing the rows indexed by / and the 
columns indexed by J from Mq, and the subscript K indicates that we are setting the 
variables a indexed by K to 0. These two generalizations are related - every Dodgson 
polynomial can be expressed as a sum of signed spanning forest polynomials (see [7j). Thus 
we can use determinant identities to derive identities for spanning forest polynomials. For 
any square matrix M, we have the classical Dodgson identity 

det(M(12, 12)) det(M) = det(M(l, 1)) det(M(2, 2)) - det(M(l, 2)) det(M(2, 1)) 

which was popularized by Dodgson through his condensation algorithm [5]. Let G be a graph 
of the form 




with two edges labelled 1 and 2, connecting three vertices V\, t>2 and V3 from top to bottom. 
The Dodgson identity gives the spanning forest polynomial identity (see section |3| 

«■■ GGrGlG»GlG»GiGi 

where for example, the graph with labels 1, 1, 2 on the vertices v±, t>2, t>3 represents $^ with 
P = {v u v 2 }U{v 3 }. 

This result can be interpreted as saying that if we transfer an extra edge from the left 
hand factor of the left hand side to the right hand factor of the left hand side, thus cutting a 
spanning tree into two in the left hand factor and joining two of the three trees together in the 
right hand factor, then we get all pairs of spanning forests with exactly two trees. However, 
it is subtle to see that the counting matches on both sides, and seems to require chains of 
edges to be transferred, along the lines of the the combinatorial proof of the Dodgson identity 
due to Zeilberger [12]. 

Equation ([I]) and its combinatorial interpretation prompted us to investigate spanning 
forest polynomial identities of the form 




Our work resulted in such an identity (Theorem [7]) which is proved in this paper. For this 
result we cannot simply interpret a classical determinantal identity; the Jacobi identity on M 
(see Corollary [9J naturally gives a cubic identity for such spanning forest polynomials, while 
the usual Dodgson identities on submatrices of M can only relate spanning forest polynomials 
whose degrees differ by at most 2. Rather, we need to combine classical identities in nontrivial 
ways. 

The paper is organized as follows: In section [2] we will set up our definitions. In section [3] 
we will define spanning forest polynomials and give their relation to the minors of M. The 
main result itself is presented and proved in section |4| Finally, in section [5] we conclude with 
a discussion of the main result, its combinatorial interpretations, and possible extensions. 

2. Graph polynomials 

Definition 1. Let G be a graph and let Mq be a matrix built as in the previous section. 
Then we define 

^ G = det(M G ) 

By the matrix-tree theorem, \&g is independent of the choice of Mq- We will call \I/g the 
graph polynomial or Kirchhoff polynomial of G. We fix a choice of M = Mq for G. 

Definition 2. Let /, J, and K be subsets of the edges of G with \I\ = \J\. Let M(J, J)k 
be the matrix obtained from M by removing the rows indexed by edges of /, the columns 
indexed by edges of J, and setting oti = for all i G K. Then we define the Dodgson 
polynomials 

= det M(I, J) K 

When G is clear it will be suppressed from the notation. Also, if K = we may suppress it 
from the notation. 

Up to sign these polynomials are independent of the choice of M (see jl]). By definition 
it is evident that = ^g- Note that if any element of K appears in / or J then it does 
not appear in M(I, J), so setting it to zero has no effect. 

Contraction and deletion of edges is natural at the level of Dodgson polynomials. 

Proposition 3. Let G be a graph and let e$ denote the i-th edge in G. Then 

^ = *G\ ei 
*G,i = *G/ ej 

Proof. The first identity follows immediately from the matrix definition of ^ and the second 
from the sum of spanning trees definition. □ 

The all-minors matrix-tree theorem [S] tells us that the monomials of any ^ G J K result 
from spanning forests of G. For our purposes it is most useful to organize these spanning 
forests with the following spanning forest polynomials. 

Definition 4. Let P = Pi U P2 U ■ ■ ■ U be a set partition of a subset of the vertices of G. 
Then we define 

F e(£F 
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where the sum runs over spanning forests F = T\ U T2 U • • • U with k component trees 
so that the vertices of Pi are in tree Tj. We note that we are allowing trees consisting of a 
single vertex. 

The relation between Dodgson polynomials and spanning forest polynomials is given by 
the following proposition. 

Proposition 5. 



(2) *ci = E ±$ & 



G\(IUJUK) 
P 

where the sum runs over set partitions P of the end points of edges of I, J, and K with the 
property that the forests corresponding to each set partition become trees in both 

G\I/(JUK) and G\J/(IUK) 

Proof. For the full details, see Proposition 12 of [7J. To sketch the argument, equation ^ 
is a direct consequence of two facts. Let E[S] be the sub matrix of E consisting of columns 
indexed by S. First, the coefficient of a given monomial m in \l/^- J is 

E[JUKUF} T ' 
-E[I U K U F] 

where F is the forest corresponding to m (that is the edges of G\(I U J U K) which do not 
contribute to m). This fact follows directly from the form of M. Second, a square matrix 
formed of columns of E has determinant ±1 if the edges corresponding to those columns are 
a spanning tree of G, and has determinant otherwise. This fact is the matrix-tree theorem 
in its most stripped down form, see for example [I] Lemma 20. □ 

3. The classical Dodgson identity 

In this section we interpret the classical Dodgson identity in terms of spanning forest 
polynomials. Consider the graph G 



det 




Apply the Dodgson determinant identity to the matrix M for G 

det(M(l, 1)) det(M(2, 2)) - det(M(l, 2)) det(M(2, 1)) = det(M) det(M(12, 12)) 
Interpreting this in terms of Dodgson polynomials gives 

and after setting the variables for edges 1 and 2 to we obtain 



For a generalization, see Corollary 10 Using the deletion-contraction relations we obtain 



^G\e l/e2 ^G\e 2 / ei - (^) 2 = *G*G\{e 1>e2 } 
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and converting to spanning forest polynomials we find that 



H t J ~ J y H J — H H 

where H is the graph with edges 1 and 2 removed. Rearranging and cancelling the squared 
term we find that 

$i a ' b ' C} $i a} ' W ' {c} = $i a ' h >'W$i a ' c >'W , ^{a,b},{c}^{a},{b,c} ^{a,c},{b}^{a},{b,c} 

ti ti ti U Jti ti ti ti 

which is just equation Q written in the spanning forest polynomial notation. See Proposition 
22 in [7] for more details. 



4. The main result 

In section [3] we gave the spanning forest polynomial version of the Dodgson identity. The 
main result of this paper is an analogous spanning forest polynomial identity for 4 marked 
vertices. Let us specialize our notation to this situation. 

Definition 6. Let Vi, t>2, t>3, and V4 be four distinct vertices of a graph G. We will write 
(ci, 02,03,04) with Ci G {1,2,3,4, — } to denote the spanning forest polynomial of the graph 
G defined by the partition of {vi : q 7^ — } with one part for each distinct integer i in 
(01,02,03,04) defined by {vt : q = £}, and no other parts. Particularly useful are the 
following special cases 

A 1 = (1,1,2,3), A 2 = (1,2,1,3), A 3 = (1,2,2,3), 
A4 = (1, 2, 3, 1), A 5 = (1, 2, 3, 2), A 6 = (1, 2, 3, 3) 

of 3 marked vertices each, the following cases 

S x = (1,1,1,2), B 2 = (1,1,2,1), B3 = (1,2,1,1), 
B 4 = (1,2,2,2), E 5 = (1,1,2,2), B 6 = (1,2,1,2), 
B 7 = (1,2, 2,1) 

of 3 marked vertices each, and finally let 

P = (1,1, 1,1). 

The Ai and Bi are the different ways in which we can partition four vertices in 3 and 2 
sets respectively. P is just for this G with four marked vertices. 
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Theorem 7. Let G be a graph with four marked vertices. Then 

(1, 1, 1, 1)(1, 2, 3, 4) = (1 - xi - x 2 )A 4 B x + x 7 A 2 B 4 + (1 - x 3 - x 2 )A 5 B 1 

+ (1 - xi - x 4 )A 6 B 1 + x 2 A 2 B 2 + (x 3 + x 2 - x 5 )A 3 B 2 
+ (1 — xi — x e )A 6 B 2 + xiA 1 B 3 + (xi - x 7 + x 4 )A 3 B 3 
+ (xi - x 8 + x 6 )A 5 B 3 + x^A x B 4 + (x 1 - x 5 + x 4 )A 3 B 5 
+ (xi - x 5 + x 6 )A 5 B 5 + x 3 A x B Q + (x 3 + x 2 - x 7 )A 3 B 6 
+ (1 - xi - x 2 + x 8 - x 6 )A 4 B 6 + (x 2 + x 7 - x 4 )A 2 B 7 
+ (1 — xi — x 7 + xs — x 6 )A e B 6 + (xi + x 5 - x 3 )A 1 B 7 
+ (1 + x b - x 3 - x 2 - x 8 )A 5 B 7 
+ (1 — xi +x 7 — x 4 — x 8 )A 6 B 7 
+ x 8 A 4 B A + x 4 A 2 B 5 + x 6 A 4 B 5 

for any x x , ...,x 8 . 

This is the generalization of the classical Dodgson identity phrased in terms of spanning 
forest polynomials. It is possible to give a graphical representation of this identity like in 
equation Q but it would take too much space. 

Outline of proof. Here we will describe the structure of the proof and the necessary calcula- 
tions will be carried out in the results which follow this outline. 

Let E(xi, x 2 , x 3 ,x 4 , x 5 , x 6 , x 7 , x 8 ) be the right hand side of We first check that E does 
not depend on the values of the Xi by checking that the coefficient of each x% in E is zero 
(Proposition 14). Now we are free to use any choice of Xi which is algebraically convenient. 

Next, using the classical Jacobi identity on an auxiliary graph with three extra edges we 
obtain an expression for 

(1,2,3,4) 2 (1,1,1,1) 



which is a linear combination of products of the form A±Aj (Lemma 15). 

Then we calculate each PAi as a linear combination of products of the form AjBk (Lemma 

2 



16). Using this calculation we obtain an expression for ((1,2,3,4)(1,1,1,1)) , which we can 
then check is the same as E(0, 0, 0, 0, 0, 0, 0, 0)E(0, 1, 0, 1, 1, 1, 1, 1) (Lemma 0. 

The proof of the theorem concludes by checking the sign. □ 

Let us pause here for a brief word on the role of the Xj. The Dodgson identities give 
a number of quadratic identities between the A\ and Bi. Consequently there cannot be 
a unique way to write (1, 1, 1, 1)(1, 2, 3, 4) as a linear combination of products AiB^. The 
Xj's describe this nonuniqueness. We can specialize to get more manageable equations, for 
example setting all Xi = and collecting terms gives 

(1,2, 3,4)(1, 1,1,1) = (1,2,3, 1)(1,-, 1,2) + (1,2,3,2)(- 1,1,2) + (1,2,3, 3)(-,-, 1,2) 

but no such specialization is canonical, so we gave the general equation in Theorem [7} 

4.1. Preliminary results. Now we can proceed with the lemmas. We will need a particular 
form of the Jacobi determinantal identity and some further Dodgson identities which follow 
from it. 
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Let M be an n x n matrix. Let / and J be subsets of {1,2, ... ,n}. Let M(I, J) be the 
matrix obtained from M by removing rows I and columns J. Similarly let M[I, J] be the 
matrix where we only keep rows I and columns J. Finally we let 

s ( J > J ) = 

Theorem 8. Let M be a nonsingular nxn matrix and let I and J be two sets in {1,2,..., n} 
with \I\ — \J\ —t. Let A = adj M and define the matrix B by bij = det(M(i, j)). Then 

det(B[I, J}) = (det M) <_1 det(M(J, J)) 

Proof. To remain self contained we will give a proof following the idea of the proof of Lemma 
28 of [I]. Let I n be the nxn identity matrix. Then 

AM = J n (det M) 

Take determinants to get 

det (A) = detM"" 1 . 

Now if the fc-th element of I is ik and the k-th element of J is jk let C be M with the jk 
column replaced by e ik , where is the z-th standard basis element of lR n . Then multiplying 
out column by column we get that AC is the matrix D whose j-th column is 

(det M)e,j if j is not in J 
Ae ik if j is j k in J 



Now notice that 
and 



detC = (-1)^' J ) det(M(/, J)) 



det D = (det M) n "* det (A[J, I}) 

= (detM) n - i det(5[J, J])(-l) s(7 ' J) . 

The second equality holds since A[J, I] can be converted to B[I, J] by multiplying each row 
and each column which had an odd index in M by — 1 and then taking a transpose; on 
determinants this changes the sign s(I, J) times. 

Finally, taking the determinant of AC = D, using the above calculations and dividing by 
(det M) n ~ l gives us the result. □ 

This formula can readily be translated into the Dodgson polynomials language. 

Corollary 9. Let G be a graph and M be its associated matrix. Let I, J and E be subsets 
of the edges, such that |/| = \ J\ = k. Then the k-level Dodgson identity is 

where I = {ii, . . . , Ik} and J = {J\, . . . , J^}. 

Proof. Use Theorem E By definition det M = ^ G and det M(I, J) = Now B[I, J] y = 

det(M(Ji, Jj)) = m 1 ^. Finally we set a e = for e G E. 

□ 



Careful book-keeping and application of the above identity yield the following corollary. 
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Corollary 10. Let M be an associated matrix for the graph G. Let E, I, J, A and B 

be ordered sets indexing edges in G, such that \AC\I\ = \B D J\ = 0, \I\ = \ J\ = k and 
\A\ = \B\ = I. Then the modified k-level Dodgson identity is 

Auh,BU.Jj\ _ . T .AUI,BUJ / lT ,A,B N k ~ ' 



(4) det i^e) 1K . j<k = nT UJ {n% 

where I = {Ii, . . . , Ik} and J = {Ji, . . . , Jk}. 

Note that when k — 2 this gives the classical Dodgson identity. 
We will use the following rearrangement of the k = 2 case. 

Proposition 11 (Brown, |4j). Let I and J be subsets of edges of G with \ J\ = |/| + 1. Let 

a, h, x be edges indices with a ^ I, b, x ^ / U J, and x < a < b. Let S = I U J U {a, b, x}. 
Then 

/ r \ -,T.Ia,J .-T.Ibx.Jx ~ Tr Iax,Jx . Tr Ib,J ~ Tr Ix,J ~ Tr Iab,Jx 

Proof. This is equation 23 from [3]; the proof proceeds by applying the k = 2 case of Q 
three times and rearranging. □ 

We only need the signs relating Dodgson polynomials to spanning forest polynomials in 
two cases, given in the next lemma. The general formula is found in Proposition 16 of [7], 
but we give here a self contained proof of the cases we need. 

Lemma 12. Fix an order and orientation of the edges of a graph G. Suppose edges 1, 2, 
and 3 have a common vertex v. Let w\, W2, and be distinct and be the other end points 
of 1, 2, and 3, and let 

if i and j are both oriented into v or both oriented out of v 
otherwise 



ijrl.2 _ £ n 2)$M>{«'i> u '2} 

V]/M = e ^^j^_iy-j+'i-^{v},{w i ,w j },{w k } 

where {i,j,k} = {1,2,3} in some order. 

Proof. The first statement of the lemma follows from the second with k = 3 applied to the 
graph G with a new vertex w 3 added and a new edge 3 from v to W3. Consider the second 
statement. Let x be the vertex which was removed when forming M. We choose it to be 
disjoint from {v,Wi,Wj}. 

Note that {v}, {wi,Wj}, {wk} is the only set partition compatible with ■ From the 
observations preceding this lemma, if ^ k ,J = then there are no common spanning trees of 
G\i/{j, k} and G\j/{i, k} and so in particular there are no terms in ${ v }<{ w i< w j}>{ w k} . Thus 

By g we know that ¥ k d = /$W.{«*.»iM«».} for some / G {-1, 1}, so it suffices to consider 
one term of Pick a term t where the tree out of uii and Wj intersects x. Let F be the 




forest corresponding to t. The sign of t in is det N where 



[ E[{i,k}UFF 
-E[{j,k}UF] 

Let B = E[{k} U F]. Then E[{i, k} U F] and E[{j, k} U F] are formed by inserting the zth 
and jth columns respectively of E into B. If {i, j} = {1,2} the insertions are both made in 
the first column. Let il be the index of the inserted column % and j' the index of the inserted 
column in j. Thus if {i,j} = {1,2} then i' = j' = 1; if {i,j} = {1,3} then {i',f} = {1,2}; 
and if {i,j} = {2, 3} then i' = f = 2. 

Consider B with the row corresponding to v removed. This is the same as the columns 
corresponding to edges of {k}UF in the incidence matrix of the graph with v and x identified. 
This has determinant ±1 since {k} U F was chosen to be a tree in this graph. Likewise, 
removing the row corresponding to W\ or w 2 we get a zero determinant since {k} U F is not 
a tree in the graph with w\ or u>2 identified with x. 

Thus if we expand det E[{i, k} U F] down the inserted column, only the cofactor coming 
from row v is retained, and likewise for E[{j, k} U F]. Thus 

det N = det(F[{z, k} U F}) det(F[{j, k} U F}) 
= e l/ e ] A-iy' +f+2e det(B) 2 

where £ is the index of row v , B is B with row v removed and e rjS is the (r, s) entry of E. □ 



4.2. Results for the main argument. Here is a catalogue of the instances of the Dodgson 
identity which we will need in the main argument, written in terms of the Aj and Bi from 
Definition |6l 

Lemma 13. 



(6) 


Ax(B 3 + B 7 ) 


+ A 2 (B 7 


-B s ) 


- M{B X + B 5 ) 


= 


(7) 


At(B 4 + B 7 ) 


+ A 5 (B 7 


-B B ) 


-A 3 (B 2 + B 5 ) 


= 


(8) 


A 2 (B 2 + B 7 ) 


+ A 1 (B 7 


-B 6 ) 


- A 4 (B X + B 6 ) 


= 


(9) 


A 2 {B 4 + B 7 ) 


+ A 6 (B 7 


-B 6 ) 


-A 3 (B 3 + B 6 ) 


= 


(10) 


A 3 (B 2 + B 6 ) 


+ A 1 (B 6 


-B 7 ) 


-A 5 (B! + B 7 ) 


= 


(11) 


A 3 (B 3 + B 5 ) 


+ A 2 (B 5 


-B 7 ) 


-A 6 (B 1 + B 7 ) 


= 


(12) 


A 4 (B 4 + B 5 ) 


+ A 5 (B 5 


-B 7 ) 


-A 6 (B 2 + B 7 ) 


= 


(13) 


A 4 (B 4 + B 6 ) 


+ A 6 (B 6 


-B 7 ) 


-A 5 (B 3 + B 7 ) 


= 


(14) 


A 5 (B 3 + B 5 ) 


+ A 4 (B 5 


-B 6 ) 


-A 6 (B 2 + B 6 ) 


= 
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Proof. The equations differ only by permuting the four marked vertices, so it suffices to prove 
(j7j). Consider the graph 




x 



1, a = 2, b = 3, I = and J = {2}, and by Lemma 12 



we 



We use identity ^ with 
obtain 

(1, -, 2, 3)(1, 2, 2, -) - (1, -, 2, -)(1, 2, 2, 3) = (1, 2, 3, 2)(1, -, 2, 2). 
For the sign of (1, — , 2, 2) note that the cutting happens first so that edges 1 and 3 become 



adjacent columns in the cut matrix. Expanding, (1,— ,2,3) 
5 4 + B 7 , (1, -, 2, -) = B 2 + B 4 + B 5 + B 7 , and (1, -, 2, 2) = 
in and rearranging gives us equation ([7]). 



= A 1 + A 3 + A 5 , (1,2,2,-) = 
£?4 + B 5 . We substitute these 

□ 



Proposition 14. yl// i/ie free variables in ^ are explained by Dodgson identities. 

Proof. The coefficient of x 3 in equation (|3| is the right hand side of equation ( [To]) , and thus 
is 0. Similarly the coefficients of £4, 25, X6, £7, and a;s are zero by fliT]), ([7]), (|14[), ([9]), and 



(13) respectively. The coefficient of £2 is in a different form, but is also zero as it is the sum 



of the right hand sides of (|8j) and (10). Finally, the coefficient of x\ is the sum of the right 
hand sides of (14), (11), and (|6|) and so is zero. □ 



Lemma 15. 



;i, 1,1, 1)(1, 2,3,4)' 




-A 3 
A 2 + A 3 + A 6 
-A* 




Proof. Let H be G with three new edges 1, 2 and 3 connecting vertex v± with the other 3 
marked vertices. By Corollary M with k = 3 and I = J = E = {1, 2, 3} we have 



(15) 



where ^> 



^ x2 T 123,123 



1,1 

H, 23 

I, 2 

H, 3 

I, 3 
H,2 



1,2 
H,3 
2,2 
i?,13 
2,3 
HA 



1,3 

H i? 
2,3 

H,l 
3,3 

H,12/ 



123,123 
H 



123,123 
P = 



is the graph polynomial of G with the edges 1,2 and 3 removed, namely 
(1, 1, 1, 1); \1/h,i23 is the spanning forest polynomial of G where each of the 



four vertices is in a separate tree, namely ^#,123 = (1,2,3,4). 

The Dodgson polynomials on the main diagonal are just spanning forest polynomials of 
G where one of the edges is removed and the other two contracted. By inspection, these are 
precisely the terms in the diagonal of the matrix in the result. The Dodgson polynomials 
on the off-diagonals require more care. We orient the edges like this: edge 2 goes towards 
vertex 1 and the other two away from it. 
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This ensures all the off-diagonal signs are negative (by Lemma 12) and that each Dodgson 
polynomial gives the desired spanning forest polynomial. The result follows. □ 



Note that the matrix in Lemma 15 is the Laplacian matrix with row and column 1 removed 
for the following graph 




where the edge labels are the A's. This is not a coincidence and there is a general identity 
which we leave out for brevity. However, the statement is analogous. 

To complete the calculation we need to multiply the whole expression by P and use the 
following 

Lemma 16. 



PA X = B X B 2 + B x B h + B 2 B 5 
PA 2 = B 1 B 3 + BtB 6 + B 3 B 6 
PA 3 = B X B 4 + B t B 7 + B±B 7 
PA4 = B 2 B 3 + B 2 B 7 + B 3 B 7 
PA 5 = B 2 B 4 + B 2 B 6 + B 4 B 6 
PA 6 = B 3 B 4 + B 3 B 5 + B 4 B 5 
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B^Bq 
B5BQ 

b$b & 

B5BQ 
B5BQ 
B$Bq 



B 5 B 7 — B e B 7 
B 5 B 7 + B 6 B 7 
B 5 B 7 + BqB 7 
B 5 B 7 + B e B 7 
B 5 B 7 + B e B 7 
B 5 B 7 — BqB 7 



Proof. By symmetry of the four vertices it suffices to prove the formula for PA\. Consider 
the graph 




Then 



PAx 



.^123,123^1,3 byLemma [l2 



= ^12,32^13,13 _ ^12,31^13,23 by Q ^ A={l},B = {3}, / = {2, 3}, 

J = {1,2}, and E = {1,2,3} 
= (1, 1, 2, -)(-, -, 1, 2) - (1, - 2, 1)(-, 1, 2, 1) by Lemma[l2] 
= (B 2 + B 5 )(B! + 52 + ^6 + 5 7 ) - (B 2 + S 7 )(B a + B 6 ) 
= B 4 B 2 + -B1-B5 + B 2 B 5 + -B5-B6 + B 5 B 7 — B e B 7 

Now we find out what happens when we multiply the equation in Lemma [15] by P. 
Lemma 17. 

((1, 1, 1, 1)(1, 2, 3, 4)) 2 = E(0, 0, 0, 0, 0, 0, 0, 0)E(0, 1, 0, 1, 1, 1, 1, 1) 
where E(x 4 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ) i/ie ng/ii /ianc? side 0/ ([3]). 
Proof. By definition 

(16) £(0, 0, 0, 0, 0, 0, 0, 0) = (A + 40 (fii + B 7 ) + A 6 (B 2 + B 6 ) + A^B, + B 6 ) 
and 

(17) E(0, 1, 0, 1, 1, 1, 1, 1) = (A + A 2 )(B 4 + B r ) + A 2 (5 2 + B 5 ) + A 4 (5 4 + B 5 ). 



□ 



Use Lemma 15 and 16 to calculate ((1, 1, 1, 1)(1,2,3, 4)) . With some trial and error we 



chose which lines of Lemma 16 to use so that the final result would look as much as possible 

\2 



like the product of (16) and (17). The term ((1, 1, 1, 1)(1, 2, 3, 4)) equals 

(At + A 2 )(PA 3 )(A 5 + A 6 ) + A 2 (PA 1 ){A 5 + A 6 ) + (A, + A 2 )(PA 5 )A G 

+ A 4 {A 4 + A 2 + A 5 + A 6 )(PA 3 ) + A 4 (A 2 + A 6 )(PA 1 + PA 5 ) 

= {At + A 2 )(A 5 + As){B x B 7 + B 4 B 7 + B X B 4 ) 

+ A 2 (A 5 + A 6 )( J B 1J B 2 + B 2 B 5 + B X B*> + B 5 B 7 ) 

+ [Ax + A 2 )A 6 (B 2 B 6 + B 4 B 6 + B 2 B 4 + B 6 B 7 ) 

+ A 4 (A 4 + A 2 + A 5 + A 6 )(B 4 B 7 + B 4 B 7 + B X B 4 + B 5 B 7 + B 6 B 7 ) 

+ A 4 (A 2 + A 6 )(B 1 B 5 + B 2 B b + B X B 2 + B 2 B d + B 4 B 6 + B 2 E 4 + B 5 B 6 ) 
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+ B 5 B 6 A 2 A 6 - B 5 BqAiA 5 - B 5 B 6 A 4 (A 1 + A 5 ) 
+ B b B 7 {A x + A 2 )A 5 + B 6 B 7 A 1 {A 5 + A 6 ) 



- B 6 B 7 
Bl) 



Now we consider the difference between this expression and (16) times (17) 
A.A^-B 2 - B 5 B 6 + B 5 B 7 + B Q B 7 ) + A 1 A 6 (-B 7 1 + B 2 B 6 - B 2 B 7 
+A 2 A b (-B 2 + B 2 B 5 - B 2 B 7 + B b B 7 ) + A 2 A 6 {-B 2 7 - 2B 2 B 7 - 
+ A 4 A 5 (B 1 B 7 + B 6 B 7 - B v B h - B 5 B 6 ) + A 4 A 6 (B 1 B 7 + B 6 B 7 + B V B 2 + B 2 B % ) 
+ A X A 4 {B 4 B 7 + B 5 B 7 - B 6 B 4 - B 5 B G ) + A 2 A 4 (B 4 B 7 + B 5 B 7 + B 2 B 5 + B 2 B 4 ) 
-Al{B x +Bs){B A + B5) 
= -{A 6 (B 2 + B 7 ) - A,(B 5 - B 7 )) {A 2 (B 2 + B 7 ) - A 1 (B 6 - B 7 )) 

- AJ^{B 1 + £? 6 )(£? 5 - B 7 ) + A 4 A & {B X + B 6 )(B 2 + B 7 ) 

- A X A 4 {B 4 + £? 5 )(5 6 - B r ) + A 2 A 4 (B 2 + B 7 )(B 5 + B 4 ) - A 2 ^ + B 6 ){B 4 + £? 5 
= -A 4 (B 4 + B 5 )(A 2 (B 2 + B r ) - A 1 (B G - B 7 )) by @ 

-jMb(£i 



£? 6 )(£? 5 - £7) + ^(Bi + £? 6 )(£? 2 + B 7 ) 
- A,A 4 (B 4 + B 5 )(B 6 - B 7 ) + A 2 A 4 (B 2 + B 7 )(B 5 + 5 4 )-^( J B 1 
A4(Bi + B 6 )(A 6 (B 2 + ,87) + A 5 (fl 7 - B 5 ) - A 4 {B 4 + B s )) 
by p| 



B 6 )(B 4 + B 5 



□ 



We are now ready to finish the proof of the main theorem. 

Proof of Theorem [?| By Lemma [17] we know that 

£7(0, 0,0, 0,0, 0,0, 0)£7(0, 1,0, 1,1, 1,1,1) = ((1, 2, 3, 4)(1, 1, 1, l)) 2 

and by Proposition 14 we know that £7 does not depend on the x^. Thus we have 

E(x u x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ) = ±(1,2,3,4)(1, 1, 1, 1) 

It remains to check the sign. Note that (1, 1, 1, 1) = P = and (1, 2, 3, 4) is \l/ for G with 
v\i v 2 , t>3, and v 4 identified. Since both (1, 1, 1, 1) and (1,2,3,4) are Kirchhoff polynomials 
of graphs all monomials appear with nonnegative coefficients. Looking at (16) we see that 
the sign is 1 and the proof is complete. 



□ 



5. Conclusions 



Theorem [7] gives a nice generalization of ([T|. Equation ([!]) itself is crucial to the combi- 
natorial and algebro-geometric approach to understanding the periods of Feynman integrals 
[TT| [2j [7J |6j [1] . In such work, having a good intuition of how to massage the polynomials which 
occur is crucial, and it is the second author's experience that spanning forest polynomials 
and their identities are very useful in this regard. 

We can ask for an edge-transferring interpretation of Theorem [7] comparable to what we 
discussed for ([I]) in the introduction. Consider (16), which is the result of setting the free 
variables to in our main theorem. Collecting terms this gives 

(1,2, 3,4)(1, 1,1,1) = (1,2,3, 1)(1,-, 1,2) + (1,2,3,2)(- 1,1,2) + (1,2,3, 3)(- - 1,2) 
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which says that we can choose to transfer an edge from any spanning forests contributing 
to (1,1,1,1) to one of those contributing to (1,2,3,4), so that we always merge the tree 
of the last vertex from (1,2,3,4) into one of the other trees, and always split the last and 
second last vertices of (1,1,1,1) into separate trees. Furthermore, the identity describes 
precisely how the split trees will interact with the remaining vertices. We know of no direct 
combinatorial proof which follows this interpretation. 

We initially obtained ^ by a numerical calculation. We first picked a graph on which 
to perform the calculations - we picked K4, K$ and K §. Then we calculated each Ai and 
Bi on this graph and then formed all possible products of A's and i?'s and formed the sum 
Y2 S t x s tA s B t , where x st is a constant, 1 < s < 6 and 1 < t < 7 for a total of 42 constants, 
and solved the linear system. The initial numerical calculation could, a priori, have had 
spurious degrees of freedom, but it could not miss any true identity of the desired form. 
Consequently, ^ is the most general quadratic formula involving 4 marked vertices. 

A natural questions is what do formulae for more marked vertices look like. Numerical 
calculations show that for 5 and 6 marked vertices the formulae have 15 and 24 free variables. 
For the classical Dodgson identity, the A's and B's are the same. If we treat the A's and the 
5's as different, we have a formula with 3 free variables. Trivially, a formula for 2 marked 
vertices has no free variables. For n = 2, 3, 4, 5 and 6 the identities so far point to expressions 
having 0, 3, 8, 15 and 24 variables in formulae for n marked vertices. These numbers are 
generated by n(n — 2) for n = 2, 3, 4, 5 and 6. 
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